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Tiny neutrino masses can arise naturally via loop diagrams. After a brief review of the radia- 
tive mass generation mechanism, I present a new model wherein TeV scale leptoquark scalars 
induce neutrino masses via two-loop diagrams. This model predicts the neutrino oscillation pa- 
rameter sin 2 0i3 to be close to the current experimental limit. The leptoquarks are accessible to 
experiments at the LHC since their masses must lie below 1.5 TeV, and their decay branching 
ratios probe neutrino oscillation parameters. Rare lepton flavor violating processes mediated by 
leptoquarks have an interesting pattern: jl — » ej may be suppressed, while jl — >> 3e and jl — e 
conversion in nuclei are within reach of the next generation experiments. Muon g — 2 receives 
new positive contributions, which can resolve the discrepancy between theory and experiment. 
New CP violating contributions to B s — B s mixing via leptoquark box diagrams are in a range that 
can explain the recently reported dimuon anomaly by the D0 collaboration. 
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The standard paradigm for explaining tiny neutrino masses is the seesaw mechanism, which 
generates an effective dimension-5 operator G\ = (LLHH) /M, suppressed by the mass scale M of 
the heavy right-handed neutrino. (L here denotes lepton doublets, while H is the Higgs doublet.) 
Oscillation data suggests that in this scenario M ~ 10 14 GeV, which is well beyond the reach of 
foreseeable experiments for direct scrutiny. An interesting alternative to the high scale seesaw 
mechanism is radiative mass generation. The smallness of neutrino masses can be understood as 
originating from loop and chirality suppression factors. The scale of new physics can naturally be 
around a TeV in this scenario. The simplest among this class of models is the Zee model [|T|] where 
neutrino masses are induced as one-loop radiative corrections arising from the exchange of charged 
scalar bosons. The effective lepton number violating operator in this model is &2 — LLLe c H/M. 
To convert this operator to neutrino mass, a loop diagram is necessary, as shown in the first diagram 
of Fig. 1. Here <I>i is a charged scalar singlet transforming as (1,1,1) of SU (3) c x SU (2) L xU(l)y 
gauge symmetry and coupling to lepton doublets as fijL{L^\ with f T = —f. <I>2(1,2, — 1/2) is a 
second Higgs doublet that generates charged lepton masses. The cubic scalar coupling <t>\<t>2H in 
the scalar potential, along with the term (&2H) 2 ensure that lepton number is explicitly broken. The 
neutrino mass in this model is given by m v ~ (fmj + m}f T )/(l6n 2 Mct>), which for = 1 TeV 
and / = 10~ 3 yields m v ~ 0.05 eV, of the right order to explain atmospheric neutrino oscillations. 
The simplest version of the Zee model is however excluded by neutrino oscillation data, since it 
predicts all the diagonal entries of the neutrino mass matrix to be zero, which is inconsistent. 

In a second class of models, neutrino masses arise as two-loop radiative corrections [Q] via 
the exchange of a singly charged scalar <E>i(l, 1, 1) and a doubly charged scalar ^2(1, 1,-2), as 
shown in the second diagram of Fig. 1. The Yukawa couplings fijLiL^\ + g i je c i e C j<£>2~ , with 
f T = — /, g T = g, along with the cubic scalar coupling <5^<!>^<I>2 ensure lepton number violation. 
The effective operator of this model is Gg — LLLe c Le c /M 2 , which requires two-loop dressing to 
convert to neutrino mass. Since m v ~ (fm£gm£f T )/[(l6n 2 ) 2 M<s>] in this model, for / ~ g ~ 0.1, 
and M<£ ~ 1 TeV, m v ^0.05 eV is generated. This model is consistent with neutrino oscillation 
data, and predicts the lightest neutrino to be nearly massless. Phenomenology of this model has 
been studied in Ref. The cross section for the production of a 1 TeV ^ at the LHC (y/s = 14 
TeV) is about 20 fb, which should be observable with its decay into same sign dileptons. 




Figure 1: Loop diagrams generating small neutrino masses in the Zee model (left) and in the model of Ref. [|2|] (right). 

A classification of low-dimensional effective AL = 2 lepton number violating operators that 
can lead to neutrino masses has been given in Ref. [Q]. The list of operators includes G3 = 
L l VQ k d c H l £ikBju which appears in the context of /^-parity violating supersymmetry. The operator 
G% = L t e c u c d c Hj£ l j is the subject for the remainder of this paper, which leads to an interesting 
neutrino mass model [Q]. G% is most directly induced by the exchange of scalar leptoquarks (LQ). 
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The order of magnitude of m v arising from is m v ~ {m t mbm T fiv) /[(\6k 2 ) 2 M^^\, where \i is 
the coefficient of a cubic scalar coupling, and v = 174 GeV is the electroweak VEV. In order to 
generate m v ^0.05 eV, it is clear that M L q must be of order TeV, which would be within reach 
of the LHC. The scalar sector consists of the leptoquark multiplets £2(3,2, 1/3) = (ft) 2 / 3 , ft) -1 / 3 ) 
and # _1 / 3 (3, 1,-2/3). Assuming global baryon number conservation, the Lagrangian relevant for 
neutrino mass is 

J^ v = YijiVitf®- 1 / 3 -btf® 2 ^ (1) 

The cubic scalar coupling will generate mixing between ft) -1 / 3 and we denote the mass 

eigenstates X a , their masses and the mixing angle G. Neutrino masses are induced via Fig. [| 



Xa 




Figure 2: Two-loop diagrams contributing to neutrino mass generation via 0%. 
Combining constraints from flavor changing processes, we find the neutrino mass matrix to be 
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Here* = F* 3 /F* 3 , y = Y 13 /Y 33 , z = Y 23 /Y 33 , w = (F* 2 /F* 3 )(Y 32 /Y 33 )(m c /m t )(m s /m b )(I jk2 /I jk3 ) and 
mo — (6^ 2 sin20F 3 * 3 }33 Ijk 3 )(m t mi,m r )/[(\6K 2 ) 2 M 2 ]. Ijt 3 denotes the two loop integral function 
shown in Fig. 3, with the internal up-type quark being the top. Since the (1,1) entry is zero, and 
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Figure 3: The integral 1^ versus M\ (left); |jc|, |v|, |z| versus 8 (middle); fi — > 3e branching ratio versus m w 2/3 (right). 

w is highly suppressed, the determinant of M v is nearly zero. This leads to the predictions m\ ~ 0, 
and tan 2 613 ~ m^jm^vn 2 0\2 in the standard parametrization of neutrino mixing. This leads to 
sin 2 0i3 = (0.044 — 0.051), which is near the current limit. A consistent fit to global oscillation 
parameters is obtained. The parameters x, y, z for such a fit are plotted as functions of the unknown 
CP violating phase 5 in Fig. 3 (middle panel). We see that |x| » 1 and \y\, \z\ ~ 1. These values 
fix the branching ratios of the leptoquarks: T(cq 2 / 3 e + b) : T(cq 2 / 3 ji + b) : T(cq 2 / 3 x + b) = 
\y\ 2 : \z\ 2 : 1, and T(X~ 1/3 -> fx~t) : r(X~ 1/3 t~t) = \x\ 2 : 1. Measuring these decays will 
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thus probe CP violation in neutrino oscillations. From the experimental limit on \i —> 3e and fi — e 
conversion in nuclei, we derive the upper limit on |F* 3 F23| as a function of ft) 2 / 3 mass. Since 
neutrino oscillation data requires Y% to be of the same order for / = 1 — 3, one can also determine 
an upper limit on F33. Combining these, we obtain an upper limit of 1.5 TeV on Mi, as shown in 
Fig. 3 as a function of \i — >► 3e branching ratio (right panel). 

Since the leptoquarks of the model are at the TeV scale, they can mediate lepton flavor vi- 
olation through the and couplings. Experimental limits on these couplings have been 
satisfied in our neutrino fit. The processes we consider are ji~ — » e~y, ji~ — » e + e~e~, \i — e 
conversion in nuclei, T~ —> e~T], T~ — >> fi ~T], B s ^ — B s ^ mixing, K — K mixing, D — D mix- 
ing, Df i^v decay, muon g — 2, 7i + jU + v^ decay, and neutrinoless double beta decay. In- 
terestingly, we find that fi ey mediated by the CO LQ is suppressed by a GIM-like mecha- 
nism. For the % LQ mass of 1 TeV, we obtain from ji — » ey, \T,i^uPzi\ < 6.7 x 10~ 3 , but no 
constraint for the CO LQ. If all the LQ masses are 1 TeV, we obtain from ji —> 3e the limits 
IF13F23I < 7.6 x 10~ 3 , IF13F23I < 1.8 x 10~ 3 , and from fi — e conversion in nuclei slightly bet- 
ter limits IF13F23I < 4.6 x 10" 3 , |Fi 3 F 23 | < 1.9 x 10" 4 . The decays z~ — >• e~v\ and t" \i~r\ 
provide the limits |Fi 2 F 32 | < 1.2 x 1(T 2 , |F 2 2F 32 | < 1.0 x 1(T 2 , now for the CO LQ mass of 300 
GeV g — 2 of the muon receives new positive contributions from the % LQ, which can be as large 
as 8(g — 2) « 12 x 10~ 10 for % mass of 300 GeV. This will be nicely consistent with the indication 
that 8(g-2) = (24.6 ±8.0) x lO" 10 . 

Recently the D0 Collaboration has reported a 3.2 sigma excess in the like-sign dimuon asym- 
metry compared to theory. A likely explanation is that there is a new source of CP violation 
in B s — B s mixing, which can arise from leptoquark box diagrams. If the B s — B s mass dif- 
ference is written as AM S = AMf M |l +h s e 2lGs \ 9 then a good fit to the data is obtained for [§] 
{h s - 0.5, a s - 120°} or {h s - 1.8, o s - 100°}. h s - 0.5 is realized in our model for |F 32 | - 1 
and ma, = 390 GeV or IF32I ~ 0.77 and = 300 GeV. The phase o s is unconstrained. This 
scenario will predict the branching ratio for B s t + t~ at the 0.25% level, compared to the stan- 
dard model value of 10~ 6 . Finally, neutrinoless double beta decay proportional to neutrino mass 
is suppressed in this model. However, it can proceed via the vector-scalar exchange process [^]. 
The diagram involves exchange of one leptoquark and one W boson. We obtain the constraint 
|F*Fn I < 1.7 x 10- 6 (^L) 2 (^P) from this process, indicating that neutrinoless double beta 
decay may be observable, in spite of the mass hierarchy being normal. 
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